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The one- and two-orbital double-exchange models for manganites are studied using Monte Carlo 
computational techniques in the presence of a robust electron-phonon coupling (but neglecting the 
antiferromagnetic exchange Jaf between the localized spins). The focus in this effort is on the 
analysis of charge transport. Our results for the one-orbital case confirm and extend previous 
recent investigations that showed the presence of robust peaks in the resistivity vs. temperature 
curves for this model. Quenched disorder substantially enhances the magnitude of the effect, while 
magnetic fields drastically reduce the resistivity. A simple picture for the origin of these results is 
presented. It is also shown that even for the case of just one electron, the resistance curves present 
metallic and insulating regions by varying the temperature, as it occurs at finite electronic density. 
Moreover, in the present study these investigations are extended to the more realistic two-orbital 
model for manganites. The transport results for this model show large peaks in the resistivity vs. 
temperature curves, located at approximately the Curie temperature, and with associated large 
magnetoresistance factors. Overall, the magnitude and shape of the effects discussed here closely 
resemble experiments for materials such as Lao.roCao.aoMnOs, and they are in qualitative agreement 
with the current predominant theoretical view that competition between a metal and an insulator, 
enhanced by quenched disorder, is crucial to understand the colossal magnetoresistance (CMR) 
phenomenon. In spite of this success, it is argued that further work is still needed to fully grasp 
the experimentally observed CMR effect, since in several other Mn oxides an antiferromagnetic 
charge-ordered orbital-ordered state is the actual competitor of the ferromagnetic metal. 

PACS numbers: 75.47.Lx, 75.30.Mb, 75.30.Kz 



I. INTRODUCTION 



One of the most outstanding open problems in the 
area of transition metal oxides is the explanation of the 
colossal magnetoresistance (CMR) effect that appears in 
the Mn oxides that are widely referred to as manganites. 
These compounds present a rich phase diagram with 
a variety of competing states which are stabilized by 
changing the carrier concentration using a standard 
chemical doping process involving ions with different va- 
lences, or by varying the carrier bandwidth via isovalent 
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Notorious among the low-temperature regimes stabilized 
in manganites are a ferromagnetic (FM) metallic phase 
and several antiferromagnetic/charge/orbital ordered 
insulating states. For the compounds with intermediate 
or small Curie temperatures, the experimentally ob- 
tained resistivity vs. temperature curves present a sharp 
peak, which occurs precisely at the transition toward 
ferromagnetism. In the vicinity of this peak, the CMR 
effect is observed, which consists of enormous changes in 
the resistivity upon the introduction of relatively small 
magnetic fields. Although technological applications of 
CMR compounds in the read-sensor industry will still 
need an increase by at least a factor two of the currently 
available critical temperatures where the large magneto- 
effects occur, the physics behind this remarkable CMR 



phenomenon defines a challenging basic-science problem 
that has attracted the attention of the condensed matter 
community. 

The explanation of the CMR effect is certainly the cru- 
cial goal of theoretical investigations in the manganite 
context. Early theoretical work showed that the stan- 
dard double exchange (DE) model was not sufficient to 
understand these materials^^. In fact, a DE model can- 
not even produce an insulator at high temperatures, in 
the realistic regimes of electronic densitiesSS, and this 
pointed toward the importance of other couplings, such 
as electron-phonon, for a proper description of these 
^iS&apoundsi^. Progress was later made with the real- 
ization that manganite models have tendencies toward 
mixed phase regimes, typically involving metallic and 
insulating states in coexistenc o"^i^° . This discovery was 
possible only after the DE model and its close variations 
were studied with unbiased computational methods be- 
yond mean-field approximations. Inhomogeneous states 
with a variety of length scales appear frequently in these 
studies and the full strength of computational techniques 
is clearly needed to fully understand this and other fam- 
ilies of complex oxides^. The theoretical discovery of 
phase separation tendenciesi^ triggered an enormous ex- 
perimental effort that confirmed the relevance of mixed 
states in most of the CMR compounds (for a review see 
Ref. 0) . Percolative pictures were envisioned to under- 
stand these materials. Model calculations by Mayr et 
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alm^ and Burgy et alMt^, using simplified spin systems 
and random resistor networks, revealed a phenomenol- 
ogy very similar to that of real CMR materials in the 
regime of couplings and electronic densities where metal- 
lic and insulating states were in competition. In fact, a 
robust peak in the resistivity of the resistor network was 
found at intermediate temperatures, and a huge change 
in its value was found to occur in the presence of magnetic 
fields. The key role of quenched disorder was remarked in 
these investigations, to obtain large enough efi'ect o'^^i'^^ . 

This initial effort using simple models was followed 
by calculations of resistivities in the more realistic, al- 
though still simplified, one-orbital model for mangan- 
ites. Verges et al2^ numerically showed that an insu- 
lator can appear at intermediate and large temperatures 
if the electron-phonon coupling A is robust enough. This 
regime, caused by localized polarons, is followed by a 
rapid transition to a metal upon cooling. In this model, 
the two tendencies in competition at low temperatures 
are both ferromagnetic, and they only differ in the char- 
acter of the charge distribution (uniform vs. localized). 
Although having two competing FM states cannot solve 
the entire CMR issue, since often the competition in 
experiments is between a ferromagnet against an anti- 
ferromagnetic/charge/orbital ordered state, the results 
were sufficiently interesting and challenging that they 
deserved further work. Moreover, they could be of rele- 
vance to important Mn oxides such as Lao.yCao.sMnOs, 
which at least naively seem well separated from charge 
ordered states in the phase diagrams. Recently, Kumar 
and Majumdar— made a very important contribution by 
proposing a new Monte Carlo (MC) algorithm to study 
fairly large lattices of the one-orbital model including 
phonons. Their main observation is that the clean limit 
results of Ref. [13 are much enhanced by including on- 
site quenched disorder, together with a robust electron- 
phonon coupling'^'^. This role of disorder to trigger po- 
laron formation in systems with strong electron-lattice 
coupling is an effect that complements the nanoscale 
phase coexistence near a first-order transition boundary 
also triggered by quenched disorder emphasized in other 
studiea^i*^. Large peaks in the resistivity vs. tempera- 
ture curves were reported in Ref. l33l resembling experi- 
ments for some manganites. Their conclusion regarding 
the importance of disorder was in agreement with previ- 
ous investigations^iSiSi^SiSi, and provided further con- 
firmation of the currently widely accepted view of man- 
ganites, namely that the essence of CMR lies in the com- 
petition of phases (metal vs. insulator), supplemented 
by quenched disorder to obtain large enough effects in 
sufficiently wide regions of parameter space. 

In spite of this tremendous progress, there are still sev- 
eral aspects of the CMR problem that need further re- 
finement. Two issues are notorious: (i) It is important 
to show that the results previously obtained for the one- 
orbital model, focusing on the resistivity peak, do also 
appear for a more realistic two-orbital model. Several 
manganites present orbital order and, as a consequence, 



using two orbitals per Mn atom is crucial for a proper 
description of these materials; (ii) The consideration of 
the antiferromagnetic spin coupling Jaf between the lo- 
calized t2g spins is also crucial. For instance, this cou- 
pling is needed to stabilize several important phases with 
charge/orbital and antiferromagnetic order, as previously 
showr.10'38. The full understanding of the CMR effect 
needs these two extra refinements. 

In the present paper, investigations of both the one- 
and two-orbital models for manganites are reported, with 
emphasis on the resistivity vs. temperature curves. The 
main results discussed in this paper are the following: (1) 
We provide an independent confirmation of the results of 
previous investigations discussed in Refs.l2ll and 33. The 
study of resistivity in Ref. relied on the analysis of 
the optical conductivity extrapolated to zero frequency. 
In addition, a novel algorithm was used to generate the 
classical spin configurations. In our present effort, a dif- 
ferent numerical method (exact diagonalization) is used 
and, more importantly, the transport properties are es- 
timated using the Landauer formalism based on trans- 
mission coefficients. Fortunately, our study shows that 
the results of Kumar and Majumdar— reporting robust 
resistivity peaks in the one-orbital model are indeed qual- 
itatively correct, even when a fairly different approach is 
used to calculate transport properties. This confirmation 
of previous investigations helps providing a robust foun- 
dation to computational studies of models for mangan- 
ites. (2) Still within the one-orbital model context, here 
a comprehensive analysis of the inffuence of quenched 
disorder and electronic density is provided, thus consid- 
erably extending the studies reported in previous efforts. 
A surprising result is that even just one Cg electron on a 
finite lattice can present transmission characteristics that 
include a resistance vs. temperature curve in qualitative 
agreement with results at finite electronic densities n and 
with experiments. Charge localization is found to be re- 
sponsible for all these features, as previously remarked 
in Rcfs. I21II33 as well. A toy example is here discussed 
to understand these results in very simple terms. (3) Fi- 
nally, one of our main new contributions is the extension 
of the previous investigations in the one-orbital model 
context into a two-orbital model framework. After a com- 
prehensive analysis of the two-orbital model properties, 
reported here, it is concluded that this model presents a 
phenomenology similar to that of the one-orbital model 
simulations, at least for the case Jaf=0. In other words, 
sharp peaks in the resistivity vs. temperature are ob- 
served in robust regions of parameter space. This con- 
clusion adds more evidence that theoretical investigations 
are on the right track toward an understanding of the 
challenging CMR effect. 

However, the important inclusion of Jaf is postponed 
for future investigations. Working at the special case 
Jaf=0 much simplifies the numerical analysis, particu- 
larly regarding the convergence properties: the localized 
spins do not have confficting tendencies, such as ferro- 
and antiferro arrangements with close energies, thus they 
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rapidly tend toward ferromagnetic ground states at low 
temperatures, even if the initial starting Monte Carlo 
configuration is random. The assumption Jaf=0 effec- 
tively reduces the global effort to merely making sure 
that the classical lattice displacements regulated by the 
electron-phonon coupling are properly converged. Impor- 
tant metastabilities were not found in our investigations. 
The considerably more subtle technical difhculties that 
will arise with the inclusion of Jaf in the model Hamil- 
tonians are left for future considerations. 

The organization of the paper is the following. In Sec- 
tion II, results for the one-orbital model are presented, 
starting with a brief discussion of the model itself and 
technical aspects. The main portion of this section is de- 
voted to the numerically calculated resistivity vs. tem- 
perature curves, obtained at several electronic densities, 
electron-phonon couplings, and strength of the quenched 
disorder. The results for the one electron problem are in- 
cluded in this section, together with evidence that charge 
localization is responsible for the insulating regime. A 
simple toy example is presented to understand the re- 
sults. In Section III, a similar analysis is presented but 
using the two-orbital model. Section IV contains the con- 
clusions of our effort, and suggestions for further work. 

II. ONE-ORBITAL MODEL 
A. Definition and Details of Simulation 

The one-orbital model used in this study is given by: 

Hlb = -t ^ (cl,QCj,a -I- h.C.) - ^ c|,Q^a,/9Ci,/3 ' Si 

2,7,0 2,7 

+ ^(A, (1) 

where cj ^ creates an electron at site i with spin a, aa,p 
are the Pauli spin matrices, (ij) indicates summing over 



where 9i and are the spherical coordinates of the core 
spin at site i (assumed classical). The operators dj now 
create an electron at site i with spin parallel to the core 
spin at i, and Ui — d\di. Note that for an infinite Hund 



nearest neighbor sites, and t is the nearest neighbor hop- 
ping amplitude for the movement of electrons (t also sets 
the energy unit, i.e. t=l in all of the results below). The 
first and second terms are the standard for a double ex- 
change model, with Si being a classical localized spin that 
represents the t2g degrees-of-freedom. The third term in 
the Hamiltonian accounts for the energy corresponding to 
the lattice-carrier interaction, with A being the strength 
of the electron-phonon coupling. Ui^^ are the distortions 
(lattice displacements) of the oxygen atoms surround- 
ing a Mn ion at site i. The index 7 in 3D (2D) runs 
over three (two) directions x,y and z {x and y). The 
tendency toward increasing the magnitude of the lattice 
distortions is balanced by the fourth term in the Hamilto- 
nian, which represents the stiffness of the Mn-0 bonds. 
Since the study of quantum phonons in this context is 
not possible with currently available algorithms, the oxy- 
gen displacements are considered classical, approxima- 
tion widely used in studies of manganites"'^''. Finally, the 
last term corresponds to the quenched disorder, which 
here it is introduced in the form of random site energies. 
Ai represents the strength of the disorder at a given site, 
and these numbers are chosen from a bimodal distribu- 
tion of width 2A with mean 0. The overall electronic 
density n is controlled with the help of a chemical po- 
tential /i added to the last term in the Hamiltonian. In 
the rest of the paper, for simplicity spatial labels will be 
denoted without arrows or bold letters independently of 
the dimension. Also the notation i -\- j is meant to rep- 
resent the lattice site given by the vectorial sum of the 
vectors corresponding to i and j, respectively. 



In this manuscript, the limit of an infinite Hund cou- 
pling will be considered, which is another widely used 
simplification known to preserve the essential physics of 
manganitesifi. In this limit, the spin of the eg-electron 
perfectly aligns along the localized t2g-spin direction, and 
the Hamiltonian is reduced to: 



(2) 

I 

coupling, the system can be shown to be particle-hole 
symmetric with respect to density n=0.5. Thus, results 
at densities n and 1-n are equivalent. 

The technique used here to handle this Hamil- 



Hlb = -t^{[cosy cos-|- +siny sin-|-e'('^'~'^^^]4dj +I1.C.} - Ai^(w, - u^^^)d\d^ 
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tonian involves the standard exact diagonalization 
of the quadratic fermionic sector for a given spin 
backgroundifi^. The procedure then consists of an evo- 
lution in Monte Carlo steps, where new spin configura- 
tions are accepted or rejected according to a standard 
Metropolis algorithm. Details have been widely dis- 
cussed in previous studies and they will not be repeated 
here^&^. Thermal averages of operators such as spin- 
spin correlation Si ■ Sj are calculated by carrying out 
an average over all Monte Carlo steps during the MC 
evolution, after discarding the initial set needed to ther- 
malize. Correlation functions at a particular distance are 
obtained by averaging over all the possible pairs of sites 
separated by that distance. As example, the definition 
of the spin correlations at distance x is the following: 



1 ■^-^ Tr[S'i • Si+xG 
Tr[e-/3^] 



(3) 

where /3 is the inverse temperature and N is the total 
number of sites, and the rest of the notation is standard. 
In the one-orbital study, mainly lattice sizes 8x8 and 
4x4x4 were used. In addition, 10'* steps were typically 
employed for thermalization, followed by another 10* for 
measurements. For larger lattices, such as 12 x 12 and 
6x6x6, 10* measurement steps were performed after 
2, 000 steps for thermalization. Most of the simulations 
have started with a random configuration of spins, but 
simulations with a FM starting configuration have also 
been carried out in order to check for convergence. No 
problems were found in this context, namely both ap- 
proaches led to very similar results. Furthermore, in- 
dependent Monte Carlo runs corresponding to different 
starting random seeds for the initial random spin con- 
figuration have also been averaged wherever possible to 
increase the accuracy of the results. 

The resistivity p has been calculated by taking the in- 
verse of the mean conductivity a, where the latter is re- 
lated to the conductance G by G = aL'^^^, with d being 
the dimension and L the linear size of the lattice. The 
calculation of the conductance G has been carried out 
following the approach extensively discussed before by 
Verges^. The use of the resistivity notation is to facil- 
itate the interpretation of results and comparison with 
experiments, namely we do not claim to have observed 
Ohmic behavior in our small system simulations. For 
the purposes of our paper, whether the resistivity or re- 
sistance is used as the key observable the conclusions are 
the same. The units used for the resistivity in the en- 
tire manuscript are [h/e^] in 2D, and [h/e'^] x L in 3D. 
Precisely in 3D, the results presented in figures were ob- 
tained by multiplying the resistance by the linear size L, 
assuming a lattice spacing one. To restore the proper 
units to our results, the real lattice spacing of Mn oxides 
must be used. 

Finally, it is important to remark that a sizable por- 
tion of the computational work presented here was car- 
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FIG. 1: (Color online) Monte Carlo results obtained using a 
4x4x4 lattice. Shown are the resistivity and spin-spin cor- 
relations, the latter at the maximum allowed distance (2\/3), 
vs. temperature, working at A=0.9, n=0.3, and for the dis- 
order strengths A indicated. The results shown are mainly 
for one configuration of quenched disorder, but as many as 10 
configurations were used in particular cases of temperatures 
and A's, and no substantial deviations were observed between 
disorder configurations. 



ried out on parallel supercomputers, in particular on the 
NCCS XT3 supercomputer (2.4-GHz AMD Opteron pro- 
cessor and 2 GB of memory) at Oak Ridge National Lab- 
oratory. Typical simulations in this effort made use of 
100 to 200 nodes in parallel. These supercomputer re- 
sources have decreased substantially the amount of real 
time that would have been needed. Indeed, we estimate 
that the entire effort would have taken at least one year to 
complete on standard small-size computer clusters. The 
message-passing interface was used to parallelize the runs 
that sweep over the various Hamiltonian parameters such 
as A and temperature. Furthermore, quenched disorder 
adds an extra level of computational effort since it re- 
quires the simulation and average of results from many 
different configurations. This extra level of complexity 
has also been parallelized. 



B. Density n=0.3 

The discussion of our computational results starts at 
the electronic density n=0.3 (equivalent to n—0.7, due to 
the symmetry discussed in the previous section). Figure 
^is a typical example of the resistivity curves obtained in 
the present effort. Shown are both the spin-spin correla- 
tion at the maximum distance possible in the cluster un- 
der study and the resistivity, working at a fixed electron- 
phonon coupling A=0.9, and varying the strength of the 
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FIG. 2; (Color online) Influence of the electron-phonon cou- 
pling A on the resistivity curves in the clean limit A=0, con- 
sidering a 4 X 4 X 4 lattice, and n=0.3. 



FIG. 3; (Color online) Effect of magnetic fields (values indi- 
cated) on the resistivity curves, for A=0.9 and A=0.6, using 
a 4 X 4 X 4 lattice. The density is n=0.3. 



quenched disorder A. In the clean hmit, A=0, there is 
a rapid change in resistivity near the transition to ferro- 
magnetism. This is a typical pure double-exchange be- 
havior: in the absence of a sufhciently strong A, quenched 
disorder, or other couplings that may lead to competing 
states, then a metal is obtained at temperatures above 
the Curie temperature. As already clearly established in 
this field, pure double-exchange models are not enough to 
address the physics of the CMR materials. However, note 
the dramatic effect of quenched disorder on the resistiv- 
ity, as shown in Fig. ^ As already recently remarked by 
Kumar and Majumdar^'^, disorder can induce a peak in 
the resistivity that much resembles experimental results 
for typical CMR materials. Even for the small systems 
studied here, the ratio of resistivities between its maxi- 
mum and minimum values is as large as ^6 for A=0.7. 
Note the correlation between the peak location and the 
temperature where ferromagnetic order appears (signaled 
in our calculations by the value of the spin-spin correla- 
tion at the largest possible distance in the cluster under 
investigation). The comparison of our results with those 
of recent publications also show that the different ways 
to study the resistivity in Ref. and here, lead to simi- 
lar qualitative data. This is an important observation, in 
view of the approximate nature of the calculations, that 
confirms the results of Ref. |33j and shows that the peak 
in the resistivity is a robust feature of the model. 

Although a variety of previous theoretical and exper- 
imental investigations have convincingly shown the im- 
portance of quenched disorder in the CMR context, nev- 
ertheless it is interesting to observe that a resistivity 
peak can also be found by varying A even in the clean 
limit A=0, as shown in Fig. 13 This observation will 
appear repeatedly in the rest of the results discussed be- 



low, namely there seems to exist a qualitative relation 
between increasing A at small A and simply increasing A 
at A=0. This fact will be exploited in the studies pre- 
sented below to simplify our task, since a simulation with 
nonzero quenched disorder needs averages over several 
disorder configurations, rendering the effort more time 
consuming than a clean-limit analysis. However, there is 
an important difference between the two cases: observ- 
ing the resistivity peak in the clean limit requires a fine 
tuning of A. For a nonzero A, the range of couplings 
with a resistivity peak is much wider (see below in the 
n—Q.l subsection for a more detailed discussion). Fine 
tuning is not compatible with the CMR effect since the 
phenomenon appears in a large number of manganese ox- 
ides, with a distribution of A's. Nevertheless, while it is 
clear that working at nonzero A and smaller than crit- 
ical A is more realistic, to the extent that the emphasis 
of a clean-limit investigation in a fine-tuned range of A is 
restricted to the vicinity of the Curie temperature, then 
both approaches appear to lead to similar conclusions. 
Finally, note that at very low temperatures the clean 
limit result shows insulating behavior, while the results 
with a nonzero disorder strength do not present such a 
feature. Although this fact establishes an interesting dif- 
ference between the two approaches, and in addition it 
is known that some manganites do present such an up- 
turn in resistivity at low temperatures^ the issue will 
not be studied in detail in the rest of the manuscript, 
since the focus of the effort is in the resistivity peak near 
the Curie temperature. The analysis of the origin of the 
low-temperature resistivity upturn in the clean limit is 
left for future work. 

The resistivity curves with peaks at intermediate tem- 
peratures (Figs.QEJl resemble the experimental data cor- 
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FIG. 4: (Color online) Influence of the disorder strength A 
(indicated) on the resistivity vs. temperature curves, working 
at A=1.2, n—0.1, and using an 8 x 8 lattice. 



FIG. 5: (Color online) Influence of the electron-phonon cou- 
pling A on the resistivity vs temperature curves, in the clean 
limit A=0, at n=0.1, and using an 8 x 8 lattice. 



responding to real manganites. It is remarkable that the 
numerical results also resemble the manner in which these 
curves are affected by magnetic fields. A typical example 
is shown in Fig. where external fields of a small value 
are used (note that these fields are small when compared 
with the natural unit, i.e. the hopping amplitude. How- 
ever, in physical units such as Teslas a field H—0.1 can 
be substantial). The region affected the most by mag- 
netic fields is the vicinity of the peak. The qualitative 
similarity with experiments is obvious, although it is fair 
to remark that in this subsection a linear scale is used for 
p while in most of the materials with truly CMR effects a 
logarithmic scale is needed, showing that the magnitude 
of the effect discussed here is substantially smaller. This 
quantitative difference could be related to the small size 
of the Monte Carlo systems used or, more fikely, with the 
absence of a strongly insulating charge-ordered orbital- 
ordered antiferromagnetic state as the direct competitor 
of the FM metallic state. Here the competitor of the FM 
metallic state is insulating but also ferromagnetic and, as 
a consequence, their resistivities near the Curie transition 
are not dramatically different. Nevertheless, even if only 
qualitatively, the similarity of the Monte Carlo data in 
Figs. Ill2i51 with experimental observations is excellent. 



C. Density n=0.1 

Several of the effects discussed at the realistic density 
n=0.3 in the previous section were found to be magni- 
fied by reducing the electronic density. In this subsection, 
the case of n=0.1 will be considered in detail (as already 
remarked, n=0.9 is the same, via symmetry considera- 



tions for the model used). The lattice to be shown is 
now two dimensional, to illustrate the similarity system- 
atically found between results in two and three dimen- 
sions. In Fig. 0] the influence of the quenched disorder 
strength A on the resistivity plots is shown. As in Fig.^ 
the case of a "small" A is considered first, namely one 
where in the clean limit the resistivity does not present 
insulating behavior. As found for n=0.3, with increas- 
ing A a prominent peak is generated, which is located at 
the Curie temperature (conclusion based on the study of 
spin correlations, not shown). The ratio of the maximum 
and minimum resistivities is now 30-50 in the range of A 
analyzed here, considerably larger than at n=0.3. 

As remarked for n=0.3, there appear to exist analo- 
gies between the processes of increasing A at "small" A 
and increasing A in the absence of quenched disorder. 
This relation is clear as well at n=0.1, and part of the 
evidence is shown in Fig. [S] which was obtained in the 
clean limit. In a narrow A range, a prominent resistivity 
peak is found, as in Fig.|21 Note the use of a logarithmic 
scale for the resistivity, showing that the magnitude of 
the effect is truly colossal. 

The effect of magnetic fields at n=0.1 is very pro- 
nounced (see Fig. E}, once again resembling the mag- 
nitude of the CMR effect in real materials. The region in 
the vicinity of the resistivity peak is the most affected. 
The magnetoresistance ratios (right panel) are as large 
as those reported in the real Mn oxides with the largest 
CMR effects. The trade-off is that the effect occurs only 
in a small window of A, but this range, as well as the mag- 
netoresistance value, can be further enlarged by adding 
quenched disorder. 

For the particular case n=0.1, it is interesting to re- 
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FIG. 6: (Color online) (a) Influence of magnetic fields on the resistivity curve and on the spin-spin correlation at the maximum 
allowed distance (4-\/2) on an 8 x 8 lattice, in the clean limit A=0. (b) Magnetoresistance ratios vs. temperature, calculated 
for two representative magnetic fields. In both (a) and (b), A=1.5 and n=0.1. 



mark the abruptness of the changes in the resistivity near 
the peak in Fig. |31 that resemble a first-order transi- 
tion The same occurs at the equivalent density n=0.9, 
as shown in Fig. [3 (a,b,c). There, the results of a longer 
Monte Carlo time simulation are presented and these 
numbers strongly suggest that indeed a first-order tran- 
sition occurs. The evidence is the jump found in (a) 
the spin correlations and (b) the resistivity. Also, in (c) 
the MC time evolution for the energy is shown. This 
presents sudden events, that resemble tunneling between 
two clearly distinct states. Note that the temperature 
chosen is slightly biased toward the highest energy state, 
since it is very difficult to fine tune T such that both com- 
peting states are visited an approximately equal amount 
of MC time. The first-order nature of the transition also 
highlights clear similarities with experiments for some 
manganites, such as Lao.vCao.sMnOa. Other Mn oxides 
appear to have a broader transition. Nevertheless, it is 
remarkable that the much simplified one-orbital model 
used here and in studies by other groups can be so rich 
to reproduce even this type of experimentally observed 
features. 

To finalize the study at density n—Q.l, it is important 
to address to what extent quenched disorder (i.e. A) does 
play a key role in generating the resistivity peak. After 
all, both in this subsection and at n=0.3 it was observed 
that even in the clean limit A=0 there is a A range where 
a peak is present. The key observation is that in the 
clean limit a fine tuning of A is needed to obtain the 
resistivity peak, namely the peak only exists in a narrow 
window of parameters. Including quenched disorder the 
range becomes much wider. For instance, in Fig. |S1 the 
area where a resistivity peak exists is shown in the A- 
A plane. This region rapidly grows with increasing A. 



Avoiding fine tuning of couplings is crucial to understand 
CMR materials, since a wide variety of Mn oxides - with 
a distribution of bandwidths and couplings - present the 
CMR effects. Any proposed mechanism must be fairly 
universal to be robust, and the inclusion of quenched 
disorder indeed renders the range of couplings for CMR 
much wider than in the clean limit. 



D. Other Electronic Densities 
and Finite-Size Effects 

The results presented thus far are only particular cases 
of the comprehensive analysis carried out in this effort, 
involving several electronic densities, couplings, and tem- 
peratures. As examples of other results obtained in the 
context of the one-orbital model, in Fig. [Tfd) results at 
n=0.7 and in the clean limit are shown. They also have 
the peak in the resistivity at the FM transition temper- 
ature, to be expected from the particle-hole symmetry 
of the model and the 71=0.3 results. A resistivity peak is 
also observed at n=0.2, as shown in Fig.jTJf), showing the 
robustness of the feature. However, the particular case 
n=0.5 is special since in this regime a staggered charge- 
ordered state is formed and, as a consequence, the system 
is strongly insulating at low temperatures in a wide range 
of electron-phonon couplings A (see Fig. [7Je)). This fact 
was also noticed in Ref. y^- Since at this density there 
is no peak in the resistivity, n=0.5 will not be further 
analyzed here. 

To complete the present analysis, size effects have also 
been investigated. It was one of our purposes to use the 
standard "exact diagonalization" (ED) method in this 
study, in order to avoid considering also issues of accuracy 
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FIG. 7: (Color online). Results mainly in the clean limit A=0 
illustrating a variety of issues discussed in the text. (a),(b), 
and (c) show the first-order (discontinuous) character of the 
transition at n—0.9, A=1.4, using an 8 x 8 lattice, (a) are 
the spin-spin correlations at the maximum distance. Also 
shown in red are results at A=0.4 and averaged over 5 disorder 
configurations, showing the smearing of the transition with 
disorder; (b) is the resistivity vs. T; and (c) is the Monte 
Carlo time evolution of the energy at T=0.034, showing the 
presence of two states, (d) is the resistivity vs. temperature 
at n—0.7 (8x8 lattice), at the A's indicated, (e) is the same as 
(d) but at n—0.5 and using a 12x12 cluster, (f) is the same 
as (e) but for n—0.2. 



if an approximate technique would have been employed, 
in addition to the intrinsic subtleties related with the 
physics involved in the problem. Moreover, we wanted 
to compare our results against those obtained with ap- 
proximate methods carried out on 8 x 8 x 8 lattices^. 
The penalization for using the ED method is that it is 
possible to carry out simulations only on up to 6 x 6 x 6 
clusters and compare those with the 4x4x4 systems 
used in the figures discussed so far. The results are in 
Fig. 1^1 The existence of the resistivity peak, the overall 
shape of the curves, and the dependence with A are very 
similar between the two lattices, supporting the conclu- 
sion that the results are robust and that indeed a CMR 
regime has been identified in these simulations (and in 
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FIG. 8: (Color online) Influence of quenched disorder on the 
size of the parameter-space region where the resistivity peak 
exists. In the plane A-A, M (I) denotes the region where the 
resistivity is metallic (insulating) at all temperatures, while 
M-I is the area where the resistivity peak is present. The 
calculation was done on an 8 x 8 cluster, with n=0.1. 



those reported before in Ref. 



E. The One-electron Problem 

The results reported in the previous sections indicate 
that the magnitude of the resistivity peak, namely the 
ratio between the maximum and minimum resistivities, 
increases when reducing the density n. In fact, the CMR 
effect is much larger at n=0.1 than at n=0.3. As a conse- 
quence, it is natural to wonder if for the case of just one 
electron a peak in the resistivity will also appear. Our ef- 
fort is carried out in the grand canonical ensemble, but it 
is possible to tune the chemical potential with sufficient 
accuracy so that just one mobile electron is MC simu- 
lated. The results are shown in Fig. ^| obtained on an 
8x8 cluster. It is remarkable to find that indeed the one- 
electron problem has a resistivity vs. temperature curve 
clearly resembling those of the other electronic densities. 
The inset of Fig. llOf a) shows that ferromagnetism in the 
classical spins is obtained in this case as well. In the bulk, 
likely only a finite-size FM region can be associated with 
a single electron, but on a finite small cluster this region 
can be as large as the entire system, as it occurs in our 
case. 

An interesting detail of the one-electron study is that 
the insulating regime is observed even at A=0. This oc- 
curs only at this very small electronic density; at n=0.1 
or 0.3, a robust value of A is needed to see a similar be- 
havior. This can be understood as follows. The cluster 
spin-spin correlations are sketched in Fig. llOf b) at low 
temperature: here the entire 8x8 cluster is ferromagnetic 
in agreement with expectations. However, at higher tem- 
peratures, in the "insulating" portion of the A=0 resistiv- 



9 




T/t T/t 

FIG. 9: (Color online) Figure to illustrate the similarity of results obtained using a 6 x 6 x 6 lattice (shown) as compared with 
the 4x4x4 results discussed before, (a) Resistivity vs. temperature at n—0.9 in the clean limit A=0, for the A's indicated, 
(b) Resistivity p vs. T at n—0.8 at the A's indicated and for A=0. 



ity curve, there are patches that are FM as well, as shown 
in Fig. llOf c). This is correlated with charge localized in 
the darker regions (not shown) .-^ Namely, in the insu- 
lating regime there is a "self-trapping" of the electrons 
that takes place, in the form of a small FM polaron. The 
lattice does not need to be distorted to see this curious ef- 
fect. At temperatures higher than r^O.125, the resistiv- 
ity now changes to a metallic state with a more uniform 
distribution of charge. The FM polarons at intermediate 
temperatures have sizes involving several lattice spacings 
and, thus, as n grows it is not surprising that their over- 
lap rapidly renders the system fully metallic. At densities 
n^O.l or larger, only with increasing A is that a charge 
localized regime (with small polarons) can be achieved. 
This issue will be discussed in more detail in the next 
subsection. 



F. Intuitive Understanding of the Results 

To intuitively understand the results, the picture 
emerging from the one electron problem is important. 
It seems that upon cooling a paramagnetic metal first 
turns into an insulator via localization of charge (this is 
the insulating regime of the resistivity curve) and then, 
fairly abruptly, a transition to a metallic FM state oc- 
curs. Although we have not calculated the entropy ex- 
plicitly (this is typically complicated to do in numerical 
simulations), we believe that in the range of A's where 
this phenomenon occurs, there is a competition between 
a FM metallic state and a charge localized (CL) state. 
The former has lower energy, but the latter has higher 
entropy due to the fact that the charge can be localized 




FIG. 10: (Color online) Results obtained in the one electron 
limit, (a) Resistivity p vs. T at A = for an 8 x 8 lattice, 
showing that even the A = results present a peak. The inset 
contains the spin-spin correlations at the maximum distance; 
(b) and (c) are the A=0 spatially resolved nearest-neighbor 
spin-spin correlations NN{i)='}2i^ij^ Si ■ Sj, where the sum is 
over the four neighbors j of site i. The results were obtained at 
r=0.004 and r=0.072, respectively, namely before and after 
the resistivity peak. Dark colors denote large values of NN{i), 
namely regions where the spins are aligned ferromagnetically. 



in a variety of arrangements. For this reason at high tem- 
perature the CL state dominates, but then a crossing to 
the FM metal occurs at low temperatures. 

This intuitive picture is compatible with a visual in- 
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FIG. 11: (Color online) (a)-(f) Density-of-states at A=1.5, 
A=0, using an 8 X 8 lattice, at the various temperatures in- 
dicated. The red lines (vertical) indicate the location of the 
chemical potential such that n=0.1 in all the panels. 



vestigation carried out in this effort. In the interesting 
coupling and density regimes, Monte Carlo snapshots of 
the classical spin configurations and electronic density of 
charge systematically reveal charge localization in the in- 
sulating region. This is correlated with the appearance 
of new structure in the density-of-states (DOS), as shown 
in Fig. ^2 In this figure, the DOS is shown for the case 
n=0.1 varying the temperature in the interesting regime 
identified in Figs. lSIOl The DOS starts developing a pseu- 
dogap (PG) feature at the chemical potential at r=0.15. 
This PG grows upon cooling and it reaches its maxi- 
mum depth at the temperature where the resistivity is 
maximized. Upon further cooling, the DOS turns into a 
typical FM density-of-states of a finite system, showing 
multiple spikesSS. The presence of a PG in the DOS of 
a model for manganites was first observed in Ref. |^ and 
our present results are compatible with those early ob- 
servations. Clearly, the resistivity peak is unrelated with 
Anderson localization that produces a mobility edge in 
the DOS, but not a PG. 

The ideas discussed here related with DOS pseudogaps 
and charge localization can be made more quantitative as 
follows. In Fig. ll2[ the inverse of the DOS at the chemical 
potential is shown, together with the logarithm of the 
resistivity. Both quantities show a similar trend with 
temperature, and the PG indeed appears correlated with 
the behavior of the resistivity. However, note that the 
use of the logarithmic scale for p indicates that the effect 
leading to the PG formation affects much more strongly 
the transport properties of the system than others. This 
is typical of a percolative system where small changes in 
the electronic distribution can lead to dramatic changes 
in the transport characteristics. 

For systems without quenched disorder the average lo- 
cal density is always constant due to translational in- 
variance. Therefore, we measure the localization of the 
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FIG. 12: (Color online) Natural logarithm of the resistivity 
p vs. temperature in the clean limit, plotted on the same 
scale with cr„ (defined in text) and the inverse of the density 
of states at the chemical potential, 1/N[uj — p), working at 
A — 1.5, n=0.1, and using an 8 x 8 lattice. Both, logj^Q(p) and 
1/N{u = /i) are normalized to coincide with the maximum 
of an- Results shown correspond to averages over several 
independent Monte Carlo runs. 



charge by calculating the error in Ui given by 



(4) 



This quantity indicates the difference between the actual 
charge rii at each site and the nominal average density 
in the full cluster, i.e. n. For a system with a uniform 
distribution of charge (t„ vanishes. This indeed occurs 
at very low temperatures. But for a system with charge 
localization then (t„ is different from zero, as it occurs 
at the resistivity peak. In Fig. ^1 cr„ is plotted versus 
temperature, showing that it follows the behavior of the 
resistivity indicating that localization of the charge is the 
main reason for the insulating regime observed above the 
Curie temperature. 

To gain further qualitative understanding for the ex- 
istence of the peak in the resistivity, we have studied a 
special cluster that seems to have common features with 
those analyzed thus far. The cluster is a 4 x 4 x 4 lattice, 
where the 32 sites on the right have a relatively small A, 
i.e. not strong enough to lead to localization of charge, 
while the 32 sites on the left have a large A. This allow 
us to clearly separate in space regions with and without 
charge localization. Monte Carlo simulating this system 
lead us to the resistivity and spin-spin correlations shown 
in Fig. I13r a') . The shape is very similar to that of other 
simulations previously described. An important point to 
notice is that the total amount of charge on the left side 
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FIG. 13: (Color online) Possible explanation for the existence of the resistivity peak, using a 4 x 4 x 4 and the couplings 
and densities indicated. Results shown are for the artificial "left-right" system described in the text, where the left (right) 
of the lattice has a relatively large (small) A. (a) Various quantities (see middle inset) vs. temperature. The growth of (T„ 
(left) and n (left) with decreasing temperature is indicative of localization of charge in the large A region (left). The trends in 
these quantities are very similar to the resistivity in its insulating range. At the Curie temperature (see spin-spin correlations), 
the localization features remains the same but now the mobile carriers can conduct much better than in a paramagnetic spin 
background, (b)-(e) Electronic density (proportional to the diameter of the spheres) at the temperatures indicated. 



of the cluster grows with decreasing temperature and the 
same does (T„, the quantity that measures the degree of 
localization. The increase of these two quantities with 
cooling (shown in the figure as well) is correlated with the 
increase of resistivity, namely with the insulating regime 
in the resistivity plot. Thus, it is clear that the insulator 
portion of the resistance is caused by charge localization. 
Examples of the local densities are in Fig. I13f b-e) . 

At the temperature where the resistivity turns metal- 
lic upon cooling, namely at the Curie temperature, note 
that the localization parameters remain approximately 
the same as at higher temperatures. Thus, the amount 
of charge that is localized does not change dramatically 
at the metal-insulator transition. What does change is 
the spin background and since the conducting properties 
of a ferromagnet and a paramagnet are very different, 
then there is a notorious reduction of the resistance be- 
low the Curie temperature. The combination of these 
two effects leads to the notorious resistivity peaks found 
in the Monte Carlo simulations. 



III. TWO-ORBITAL MODEL 
A. Definition 

The two orbitals used in this model arise from the two 
Cg bands that are active at the Mn ions in Mn-oxides, as 
extensively discussed beforeiiSiSSiS^ The Hamiltonian for 
this model i^iSiS&S^ 



H2b — ^ t'^'y'S{9i,4>i,9i+a,<t>i+a)<^i^^Ci+art' 
iPi + Q2iTxi + QsiTzi) 

i 

a=3 
i a—1 

where the factor that renormalizes the hopping in the 
Jh=oo Hmit is 

6- 9 
5(6',,0i,6'j,0j) = cos(-^)cos(-^) 

+ sin(^)sin(^)e-*('>"-^^:». (6) 

The parameters f^^, are the hopping amplitudes between 
the orbitals 7 and 7' in the direction a. In this sec- 
tion, we restrict ourselves to two dimensions, such that 
tla = -VStl, = -^/3^L = Mt, = 1, and = VStl, = 
VSt^^ = 3i^^ = 1. Qii, Q2i and Q^i are normal modes of 
vibration that can be expressed in terms of the oxygen 
coordinate u ^ as ■ 

Qli = —f^[{Ui,z - Ui-z.z) + {Ui,x - Ui^x,x) 
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Qzi — Y^^^^'^ ^i~z,z^ ^i—x,x) 
Also, Txi = ^ia^'i-b ^ '^ib'^i'^^ = cl^Cia c|^^Ci6, and 

Pi — cl^Cia + cli^Cib- The constant A is the electron- 
phonon couphng related to the Jahn- Teller distortion 
of the MnOe octahedroniSi2iiLiiiiiSLaSiSi Regarding the 
phononic stiffness, and in units of t^^ — 1, the Da pa- 
rameters are Di = 1 and D2 = D3 = 0.5, as discussed 
in previous literature4ii The rest of the notation is stan- 
dard. Note that in the large Hund coupling limit there is 
no spin index. The Jaf coupling between the localized 
spins is neglected, as for the one-orbital model. In some 
of the results below, a Zeeman term with field strength 
H was added. 

The main purpose of the numerical study discussed 
in this section is to investigate if the two-orbital model 
for manganites can also produce a resistivity peak, as ob- 
served in the one-orbital case. The study in this section is 
presented with the same caveats as the one-orbital inves- 
tigation: (i) it is an important step toward a realistic the- 
oretical description of manganites since Mn-oxides have 
two active Cg orbitals, but (ii) the model does not include 
the coupling Jaf which is crucial to generate the realis- 
tic insulating state, with antiferromagnetic and orbital 
order. Nevertheless, the observation of features that in 
several ways resemble experiments is exciting and at least 
part of the essence of real materials appears to have been 
captured by the models discussed here, even in purely FM 
regimes. We are aware that conclusions similar to ours 
have also been reached recently independently by Kampf 
and Kumar—. 



B. Density n=0.3 




T/t 

FIG. 14: (Color online) (a) Resistivity vs. T for various val- 
ues of A, using the two-orbital model for manganites. The 
simulation was performed on an 8 x 8 lattice with 20 electrons 
(n ~ 0.3), using 5,000 MC steps for thermalization and 5,000 
for measurements. These results are in the clean limit A=0. 
A ferromagnetic arrangement was used as the starting config- 
uration, although tests using paramagnetic states reveal very 
similar results, (b) Spin-spin correlation (at the maximum 
distance 4\/2 allowed in the studied lattice) vs. T for various 
values of A. Lattice, density, and MC steps are as in (a). 



Typical computational results for the two-orbital 
model at n'--^0.3 are shown in Fig. In the clean limit 
A=0, there is a narrow region of A where a well-defined 
peak is found in the resistivity. The location of the peak 
is correlated with the appearance of ferromagnetic order, 
as also shown in the same figure. The systematic ten- 
dencies and behaviors observed in the two-orbital model 
simulations are very similar to those found for the one- 
orbital model. 

Finite-size effects do not seem to modify strongly our 
conclusions, as also found for the one-orbital case. In 
Fig. El results obtained using a 12 x 12 cluster are re- 
ported. A resistivity peak is observed in a very similar 
range of A as found using the 8x8 cluster. Although these 
results are not sufficient to fully prove that the behavior 
found on finite lattice survives the bulk limit, they are 
very suggestive: in two and three dimensional lattices, 
for a wide range of electronic densities, for a variety of 
lattice sizes, with and without quenched disorder, and 
both for the one- and two-orbital models the resistivity 



peak is present in the study of charge transport. 

This resistivity peak in the two-orbital model is also 
drastically affected by relatively small magnetic fields, as 
observed for the one-orbital case. Typical results are in 
Fig. El The region the most affected by the magnetic 
fields is where the maximum in the resistivity is observed, 
as expected. The magnitude of the magnetoresistance 
effect shown in the figure is comparable to the numbers 
found for the one-orbital case at similar electronic densi- 
ties (see Fig. O. 

It is also important to discuss the influence of quenched 
disorder. Typical results are in Fig. 1171 where the resis- 
tivity vs. T is shown both with and without disorder. 
As anticipated from recent previous investigations'^'^, and 
from the one-orbital study in this manuscript, it was ob- 
served that introducing quenched disorder enhances sub- 
stantially the features found in the clean limit. For in- 
stance, at A=1.66, the A=0 curve does not show a resis- 
tivity peak, but this feature is generated at A=0.2 and 
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FIG. 15: (Color online) (a) Resistivity vs. T for a 12x12 
lattice in the clean limit A=0, considering 2 orbitals per site, 
44 electrons (n~l/3), Jaf=0.0, and the values of A indicated. 
The simulation was carried out starting with a random state 
at each temperature, warming up for 7,500 MC steps per site, 
followed by 7,500 MC steps for measurements, (b) Spin-spin 
correlation at the largest possible distance (6\/2) on a 12x12 
lattice, in the clean limit, with the same convention and pa- 
rameters as in (a). 
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FIG. 16: (Color online) (a) Resistivity vs. T for several A's 
(indicated). Results were obtained with (filled symbols) and 
without (open symbols) a magnetic field _ff=0.1. The sim- 
ulation was performed on an 8x8 lattice, with 20 electrons 
(n « 0.3), 5,000 thermahzation MC steps followed by 5,000 
MC steps for measurements, (b) Magneto-resistance (defined 
as ((p(0)-p(Ji'))/p(//))xl00 vs. T for the same parameters as 
in (a). 



the same A. In cases where the resistivity already has a 
peak in the clean limit, this structure is enlarged with 
increasing A. 



C. Density n=0.1 

In Fig.^1 Monte Carlo results at n=0.1 are presented 
for the two-orbital model. The behavior of the resis- 
tivity is very similar to what was observed for the one- 
orbital case (see Figs. I^Ja) andEJ^a)), namely a clearly 
defined peak is observed, and a sharp (likely first order) 
transition in the resistance occurs upon cooling through 
the Curie temperature. At this electronic density, the 
changes in resistance upon heating or cooling are much 
larger than at other densities such as n=0.3. 

Overall, it is clear that the models with one and two 
orbitals behave fairly similarly, and the existence of a 
peak in the resistivity is a robust result of this effort and 
previous Monte Carlo simulations^. 



IV. CONCLUSIONS 

The research effort discussed in this paper reached sev- 
eral goals. First, it confirmed recent reports by other 
groups^ii^?: regarding the existence of a large peak in the 
resistivity vs. temperature for the one-orbital model for 
manganites, including a robust electron-phonon coupling. 
This confirmation is interesting since the results of the 
previouSiS^ and current efforts were obtained using dif- 
ferent techniques to estimate transport properties, and 
also with different methods to simulate the one-orbital 
model. Second, a comprehensive analysis of the influence 
of couplings, quenched disorder strength, and electronic 
density was here described. This includes the case of just 
one electron on an otherwise carrier empty lattice, in the 
presence of classical t2g spins. This one-electron problem 
also presents a large resistance peak when varying the 
temperature. A very simple explanation for the behavior 
of these systems was discussed, based on a competition 
between tendencies to charge localization and ferromag- 
netism. 
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FIG. 17: (Color online) Resistivity vs. T for various values 
of A and with and without quenched disorder, as indicated. 
The figure shows the enhancement of the resistivity peak with 
increasing A. The lattice size, MC steps, and various param- 
eters are as in Fig. 1161 
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FIG. 18: (Color online) Results for the two-orbital model at 
n—Q.l and the A's indicated, (a) Resistivity vs. temperature 
using an 8 X 8 lattice, 5,000 MC steps for thermalization and 
a similar number of measurements. Note how sharp is the 
low temperature transition from low to high resistance, (b) 
Spin-spin correlation at the maximum distance. 



Finally, the present effort also includes a study of the 
two-orbital model for manganites. The overall conclusion 
is that its behavior is similar to that of the one-orbital 
model. Since these results are themselves also similar to 
experiments, our effort and those of other groups pro- 
vide evidence that the theoretical studies that focus on 
the regime of competition between a metal and an insu- 
lator are on the right track toward a full explanation of 
the CMR phenomenon. Both with one and two orbitals, 
quenched disorder is important to enlarge the magnitude 
of the effects and broaden its range in parameter space, 
thus avoiding the fine tuning of couplings needed in the 
clean limit. 

In the present and related investigationsSiiS both the 
metal and the insulator at low temperatures are ferro- 
magnetic, and they differ only in the arrangement of 
charge (extended vs. localized character). The next level 
of sophistication of the simulations of manganite mod- 
els must address the competition between different mag- 
netic orders. Although the current results have a remark- 
able resemblance with several experiments, it is known 
that the largest effects in real manganites occur when a 
antiferromagnetic/charge/orbital ordered state competes 
with the ferromagnetic metal. To achieve this final goal, 
the coupling Jaf must be incorporated in the investi- 
gations. This will require levels of numerical accuracy 
higher than in the present effort, due to the competition 
of very different states that typically lead to metastabili- 
ties and long thermalization times. Results will hopefully 
be presented in the near future. 
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